3.337 \(\int \frac {x^3 (c+d x+e x^2)}{a+b x^3} \, dx\)

Optimal. Leaf size=205 \[ \frac {\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )}{6 b^{4/3}}-\frac {\sqrt [3]{a} \left (\sqrt [3]{b} c-\sqrt [3]{a} d\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 b^{5/3}}+\frac {\sqrt [3]{a} \left (\sqrt [3]{a} d+\sqrt [3]{b} c\right ) \tan ^{-1}\left (\frac {\sqrt [3]{a}-2 \sqrt [3]{b} x}{\sqrt {3} \sqrt [3]{a}}\right )}{\sqrt {3} b^{5/3}}-\frac {a e \log \left (a+b x^3\right )}{3 b^2}+\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b} \]

[Out]

c*x/b+1/2*d*x^2/b+1/3*e*x^3/b-1/3*a^(1/3)*(b^(1/3)*c-a^(1/3)*d)*ln(a^(1/3)+b^(1/3)*x)/b^(5/3)+1/6*a^(1/3)*(c-a
^(1/3)*d/b^(1/3))*ln(a^(2/3)-a^(1/3)*b^(1/3)*x+b^(2/3)*x^2)/b^(4/3)-1/3*a*e*ln(b*x^3+a)/b^2+1/3*a^(1/3)*(b^(1/
3)*c+a^(1/3)*d)*arctan(1/3*(a^(1/3)-2*b^(1/3)*x)/a^(1/3)*3^(1/2))/b^(5/3)*3^(1/2)

________________________________________________________________________________________

Rubi [A]  time = 0.26, antiderivative size = 205, normalized size of antiderivative = 1.00, number of steps used = 10, number of rules used = 9, integrand size = 23, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.391, Rules used = {1887, 1871, 1860, 31, 634, 617, 204, 628, 260} \[ \frac {\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )}{6 b^{4/3}}-\frac {\sqrt [3]{a} \left (\sqrt [3]{b} c-\sqrt [3]{a} d\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 b^{5/3}}+\frac {\sqrt [3]{a} \left (\sqrt [3]{a} d+\sqrt [3]{b} c\right ) \tan ^{-1}\left (\frac {\sqrt [3]{a}-2 \sqrt [3]{b} x}{\sqrt {3} \sqrt [3]{a}}\right )}{\sqrt {3} b^{5/3}}-\frac {a e \log \left (a+b x^3\right )}{3 b^2}+\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b} \]

Antiderivative was successfully verified.

[In]

Int[(x^3*(c + d*x + e*x^2))/(a + b*x^3),x]

[Out]

(c*x)/b + (d*x^2)/(2*b) + (e*x^3)/(3*b) + (a^(1/3)*(b^(1/3)*c + a^(1/3)*d)*ArcTan[(a^(1/3) - 2*b^(1/3)*x)/(Sqr
t[3]*a^(1/3))])/(Sqrt[3]*b^(5/3)) - (a^(1/3)*(b^(1/3)*c - a^(1/3)*d)*Log[a^(1/3) + b^(1/3)*x])/(3*b^(5/3)) + (
a^(1/3)*(c - (a^(1/3)*d)/b^(1/3))*Log[a^(2/3) - a^(1/3)*b^(1/3)*x + b^(2/3)*x^2])/(6*b^(4/3)) - (a*e*Log[a + b
*x^3])/(3*b^2)

Rule 31

Int[((a_) + (b_.)*(x_))^(-1), x_Symbol] :> Simp[Log[RemoveContent[a + b*x, x]]/b, x] /; FreeQ[{a, b}, x]

Rule 204

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> -Simp[ArcTan[(Rt[-b, 2]*x)/Rt[-a, 2]]/(Rt[-a, 2]*Rt[-b, 2]), x] /
; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])

Rule 260

Int[(x_)^(m_.)/((a_) + (b_.)*(x_)^(n_)), x_Symbol] :> Simp[Log[RemoveContent[a + b*x^n, x]]/(b*n), x] /; FreeQ
[{a, b, m, n}, x] && EqQ[m, n - 1]

Rule 617

Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*Simplify[(a*c)/b^2]}, Dist[-2/b, Sub
st[Int[1/(q - x^2), x], x, 1 + (2*c*x)/b], x] /; RationalQ[q] && (EqQ[q^2, 1] ||  !RationalQ[b^2 - 4*a*c])] /;
FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 628

Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[(d*Log[RemoveContent[a + b*x +
c*x^2, x]])/b, x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]

Rule 634

Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Dist[(2*c*d - b*e)/(2*c), Int[1/(a +
b*x + c*x^2), x], x] + Dist[e/(2*c), Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x] &
& NeQ[2*c*d - b*e, 0] && NeQ[b^2 - 4*a*c, 0] &&  !NiceSqrtQ[b^2 - 4*a*c]

Rule 1860

Int[((A_) + (B_.)*(x_))/((a_) + (b_.)*(x_)^3), x_Symbol] :> With[{r = Numerator[Rt[a/b, 3]], s = Denominator[R
t[a/b, 3]]}, -Dist[(r*(B*r - A*s))/(3*a*s), Int[1/(r + s*x), x], x] + Dist[r/(3*a*s), Int[(r*(B*r + 2*A*s) + s
*(B*r - A*s)*x)/(r^2 - r*s*x + s^2*x^2), x], x]] /; FreeQ[{a, b, A, B}, x] && NeQ[a*B^3 - b*A^3, 0] && PosQ[a/
b]

Rule 1871

Int[(P2_)/((a_) + (b_.)*(x_)^3), x_Symbol] :> With[{A = Coeff[P2, x, 0], B = Coeff[P2, x, 1], C = Coeff[P2, x,
2]}, Int[(A + B*x)/(a + b*x^3), x] + Dist[C, Int[x^2/(a + b*x^3), x], x] /; EqQ[a*B^3 - b*A^3, 0] ||  !Ration
alQ[a/b]] /; FreeQ[{a, b}, x] && PolyQ[P2, x, 2]

Rule 1887

Int[(Pq_)/((a_) + (b_.)*(x_)^(n_)), x_Symbol] :> Int[ExpandIntegrand[Pq/(a + b*x^n), x], x] /; FreeQ[{a, b}, x
] && PolyQ[Pq, x] && IntegerQ[n]

Rubi steps

\begin {align*} \int \frac {x^3 \left (c+d x+e x^2\right )}{a+b x^3} \, dx &=\int \left (\frac {c}{b}+\frac {d x}{b}+\frac {e x^2}{b}-\frac {a c+a d x+a e x^2}{b \left (a+b x^3\right )}\right ) \, dx\\ &=\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b}-\frac {\int \frac {a c+a d x+a e x^2}{a+b x^3} \, dx}{b}\\ &=\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b}-\frac {\int \frac {a c+a d x}{a+b x^3} \, dx}{b}-\frac {(a e) \int \frac {x^2}{a+b x^3} \, dx}{b}\\ &=\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b}-\frac {a e \log \left (a+b x^3\right )}{3 b^2}-\frac {\int \frac {\sqrt [3]{a} \left (2 a \sqrt [3]{b} c+a^{4/3} d\right )+\sqrt [3]{b} \left (-a \sqrt [3]{b} c+a^{4/3} d\right ) x}{a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2} \, dx}{3 a^{2/3} b^{4/3}}-\frac {\left (\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {1}{\sqrt [3]{a}+\sqrt [3]{b} x} \, dx}{3 b}\\ &=\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b}-\frac {\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 b^{4/3}}-\frac {a e \log \left (a+b x^3\right )}{3 b^2}-\frac {\left (a^{2/3} \left (\sqrt [3]{b} c+\sqrt [3]{a} d\right )\right ) \int \frac {1}{a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2} \, dx}{2 b^{4/3}}+\frac {\left (\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {-\sqrt [3]{a} \sqrt [3]{b}+2 b^{2/3} x}{a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2} \, dx}{6 b^{4/3}}\\ &=\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b}-\frac {\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 b^{4/3}}+\frac {\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )}{6 b^{4/3}}-\frac {a e \log \left (a+b x^3\right )}{3 b^2}-\frac {\left (\sqrt [3]{a} \left (\sqrt [3]{b} c+\sqrt [3]{a} d\right )\right ) \operatorname {Subst}\left (\int \frac {1}{-3-x^2} \, dx,x,1-\frac {2 \sqrt [3]{b} x}{\sqrt [3]{a}}\right )}{b^{5/3}}\\ &=\frac {c x}{b}+\frac {d x^2}{2 b}+\frac {e x^3}{3 b}+\frac {\sqrt [3]{a} \left (\sqrt [3]{b} c+\sqrt [3]{a} d\right ) \tan ^{-1}\left (\frac {\sqrt [3]{a}-2 \sqrt [3]{b} x}{\sqrt {3} \sqrt [3]{a}}\right )}{\sqrt {3} b^{5/3}}-\frac {\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 b^{4/3}}+\frac {\sqrt [3]{a} \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )}{6 b^{4/3}}-\frac {a e \log \left (a+b x^3\right )}{3 b^2}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 0.11, size = 191, normalized size = 0.93 \[ \frac {\sqrt [3]{b} \left (\sqrt [3]{a} \sqrt [3]{b} c-a^{2/3} d\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )+2 \sqrt [3]{b} \left (a^{2/3} d-\sqrt [3]{a} \sqrt [3]{b} c\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )+2 \sqrt {3} \sqrt [3]{a} \sqrt [3]{b} \left (\sqrt [3]{a} d+\sqrt [3]{b} c\right ) \tan ^{-1}\left (\frac {1-\frac {2 \sqrt [3]{b} x}{\sqrt [3]{a}}}{\sqrt {3}}\right )-2 a e \log \left (a+b x^3\right )+6 b c x+3 b d x^2+2 b e x^3}{6 b^2} \]

Antiderivative was successfully verified.

[In]

Integrate[(x^3*(c + d*x + e*x^2))/(a + b*x^3),x]

[Out]

(6*b*c*x + 3*b*d*x^2 + 2*b*e*x^3 + 2*Sqrt[3]*a^(1/3)*b^(1/3)*(b^(1/3)*c + a^(1/3)*d)*ArcTan[(1 - (2*b^(1/3)*x)
/a^(1/3))/Sqrt[3]] + 2*b^(1/3)*(-(a^(1/3)*b^(1/3)*c) + a^(2/3)*d)*Log[a^(1/3) + b^(1/3)*x] + b^(1/3)*(a^(1/3)*
b^(1/3)*c - a^(2/3)*d)*Log[a^(2/3) - a^(1/3)*b^(1/3)*x + b^(2/3)*x^2] - 2*a*e*Log[a + b*x^3])/(6*b^2)

________________________________________________________________________________________

fricas [C]  time = 2.83, size = 4798, normalized size = 23.40 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3*(e*x^2+d*x+c)/(b*x^3+a),x, algorithm="fricas")

[Out]

1/36*(12*b*e*x^3 + 18*b*d*x^2 - 2*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6
+ 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)
*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e
^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*b^2*log(1/36*((-I*sqr
t(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*
c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1
/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 -
(d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)^2*b^4*d + 2*a*b*c*d^2 - a*b*c^2*e + a^2*d*e^2 + 1/6*(b^3*c^2 -
2*a*b^2*d*e)*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^
3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) +
9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a
*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2) + (b^2*c^3 + a*b*d^3)*x) + 36*b*c*x + (((-
I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*
(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1
)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e
^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*b^2 + 3*sqrt(1/3)*b^2*sqrt(-(((-I*sqrt(3) + 1)*(a^2*e^2/b^
4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^
6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/5
4*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b
)/b^6)^(1/3) + 6*a*e/b^2)^2*b^4 - 12*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/
b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d
*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^
2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*a*b^2*e + 144*a*b*
c*d + 36*a^2*e^2)/b^4) - 18*a*e)*log(-1/36*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^
3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 -
3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*
d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)^2*b^4*d - 2*
a*b*c*d^2 + a*b*c^2*e - a^2*d*e^2 - 1/6*(b^3*c^2 - 2*a*b^2*d*e)*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^
2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c
^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)
*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*
a*e/b^2) + 2*(b^2*c^3 + a*b*d^3)*x + 1/12*sqrt(1/3)*(((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)
/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^
3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/
18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*b^
4*d - 6*b^3*c^2 - 6*a*b^2*d*e)*sqrt(-(((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3
/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*
d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a
^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)^2*b^4 - 12*((-I*s
qrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*
b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(
-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3
- (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*a*b^2*e + 144*a*b*c*d + 36*a^2*e^2)/b^4)) + (((-I*sqrt(3) + 1
)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^
2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*
e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3
*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*b^2 - 3*sqrt(1/3)*b^2*sqrt(-(((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d
+ a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b
^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*
d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3)
+ 6*a*e/b^2)^2*b^4 - 12*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b
*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^
6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^
6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*a*b^2*e + 144*a*b*c*d + 36*a^2*
e^2)/b^4) - 18*a*e)*log(-1/36*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1
/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2
*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*
a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)^2*b^4*d - 2*a*b*c*d^2 + a
*b*c^2*e - a^2*d*e^2 - 1/6*(b^3*c^2 - 2*a*b^2*d*e)*((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(
-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3
- (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18
*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2) + 2*
(b^2*c^3 + a*b*d^3)*x - 1/12*sqrt(1/3)*(((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e
^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*
c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d +
a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*b^4*d - 6*b^3*c
^2 - 6*a*b^2*d*e)*sqrt(-(((-I*sqrt(3) + 1)*(a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(
b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b
^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b
^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)^2*b^4 - 12*((-I*sqrt(3) + 1)*(
a^2*e^2/b^4 - (a*b*c*d + a^2*e^2)/b^4)/(-1/27*a^3*e^3/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e
^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*d*e)*a^2*b)/b^6)^(1/3) + 9*(I*sqrt(3) + 1)*(-1/27*a^3*e^3
/b^6 + 1/54*(b*c^3 + a*d^3)*a/b^5 + 1/18*(a*b*c*d + a^2*e^2)*a*e/b^6 - 1/54*(a*b^2*c^3 + a^3*e^3 - (d^3 - 3*c*
d*e)*a^2*b)/b^6)^(1/3) + 6*a*e/b^2)*a*b^2*e + 144*a*b*c*d + 36*a^2*e^2)/b^4)))/b^2

________________________________________________________________________________________

giac [A]  time = 0.18, size = 208, normalized size = 1.01 \[ -\frac {a e \log \left ({\left | b x^{3} + a \right |}\right )}{3 \, b^{2}} - \frac {\sqrt {3} {\left (\left (-a b^{2}\right )^{\frac {1}{3}} b c - \left (-a b^{2}\right )^{\frac {2}{3}} d\right )} \arctan \left (\frac {\sqrt {3} {\left (2 \, x + \left (-\frac {a}{b}\right )^{\frac {1}{3}}\right )}}{3 \, \left (-\frac {a}{b}\right )^{\frac {1}{3}}}\right )}{3 \, b^{3}} - \frac {{\left (\left (-a b^{2}\right )^{\frac {1}{3}} b c + \left (-a b^{2}\right )^{\frac {2}{3}} d\right )} \log \left (x^{2} + x \left (-\frac {a}{b}\right )^{\frac {1}{3}} + \left (-\frac {a}{b}\right )^{\frac {2}{3}}\right )}{6 \, b^{3}} + \frac {2 \, b^{2} x^{3} e + 3 \, b^{2} d x^{2} + 6 \, b^{2} c x}{6 \, b^{3}} + \frac {{\left (a b^{6} d \left (-\frac {a}{b}\right )^{\frac {1}{3}} + a b^{6} c\right )} \left (-\frac {a}{b}\right )^{\frac {1}{3}} \log \left ({\left | x - \left (-\frac {a}{b}\right )^{\frac {1}{3}} \right |}\right )}{3 \, a b^{7}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3*(e*x^2+d*x+c)/(b*x^3+a),x, algorithm="giac")

[Out]

-1/3*a*e*log(abs(b*x^3 + a))/b^2 - 1/3*sqrt(3)*((-a*b^2)^(1/3)*b*c - (-a*b^2)^(2/3)*d)*arctan(1/3*sqrt(3)*(2*x
+ (-a/b)^(1/3))/(-a/b)^(1/3))/b^3 - 1/6*((-a*b^2)^(1/3)*b*c + (-a*b^2)^(2/3)*d)*log(x^2 + x*(-a/b)^(1/3) + (-
a/b)^(2/3))/b^3 + 1/6*(2*b^2*x^3*e + 3*b^2*d*x^2 + 6*b^2*c*x)/b^3 + 1/3*(a*b^6*d*(-a/b)^(1/3) + a*b^6*c)*(-a/b
)^(1/3)*log(abs(x - (-a/b)^(1/3)))/(a*b^7)

________________________________________________________________________________________

maple [A]  time = 0.04, size = 231, normalized size = 1.13 \[ \frac {e \,x^{3}}{3 b}+\frac {d \,x^{2}}{2 b}-\frac {\sqrt {3}\, a c \arctan \left (\frac {\sqrt {3}\, \left (\frac {2 x}{\left (\frac {a}{b}\right )^{\frac {1}{3}}}-1\right )}{3}\right )}{3 \left (\frac {a}{b}\right )^{\frac {2}{3}} b^{2}}-\frac {a c \ln \left (x +\left (\frac {a}{b}\right )^{\frac {1}{3}}\right )}{3 \left (\frac {a}{b}\right )^{\frac {2}{3}} b^{2}}+\frac {a c \ln \left (x^{2}-\left (\frac {a}{b}\right )^{\frac {1}{3}} x +\left (\frac {a}{b}\right )^{\frac {2}{3}}\right )}{6 \left (\frac {a}{b}\right )^{\frac {2}{3}} b^{2}}-\frac {\sqrt {3}\, a d \arctan \left (\frac {\sqrt {3}\, \left (\frac {2 x}{\left (\frac {a}{b}\right )^{\frac {1}{3}}}-1\right )}{3}\right )}{3 \left (\frac {a}{b}\right )^{\frac {1}{3}} b^{2}}+\frac {a d \ln \left (x +\left (\frac {a}{b}\right )^{\frac {1}{3}}\right )}{3 \left (\frac {a}{b}\right )^{\frac {1}{3}} b^{2}}-\frac {a d \ln \left (x^{2}-\left (\frac {a}{b}\right )^{\frac {1}{3}} x +\left (\frac {a}{b}\right )^{\frac {2}{3}}\right )}{6 \left (\frac {a}{b}\right )^{\frac {1}{3}} b^{2}}-\frac {a e \ln \left (b \,x^{3}+a \right )}{3 b^{2}}+\frac {c x}{b} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x^3*(e*x^2+d*x+c)/(b*x^3+a),x)

[Out]

1/3/b*e*x^3+1/2/b*d*x^2+1/b*c*x-1/3/(a/b)^(2/3)*a/b^2*c*ln(x+(a/b)^(1/3))+1/6*a/b^2*c/(a/b)^(2/3)*ln(x^2-(a/b)
^(1/3)*x+(a/b)^(2/3))-1/3*a/b^2*c/(a/b)^(2/3)*3^(1/2)*arctan(1/3*3^(1/2)*(2/(a/b)^(1/3)*x-1))+1/3*a/b^2*d/(a/b
)^(1/3)*ln(x+(a/b)^(1/3))-1/6*a/b^2*d/(a/b)^(1/3)*ln(x^2-(a/b)^(1/3)*x+(a/b)^(2/3))-1/3*a/b^2*d*3^(1/2)/(a/b)^
(1/3)*arctan(1/3*3^(1/2)*(2/(a/b)^(1/3)*x-1))-1/3*a/b^2*e*ln(b*x^3+a)

________________________________________________________________________________________

maxima [A]  time = 2.94, size = 190, normalized size = 0.93 \[ \frac {2 \, e x^{3} + 3 \, d x^{2} + 6 \, c x}{6 \, b} - \frac {\sqrt {3} {\left (a b d \left (\frac {a}{b}\right )^{\frac {2}{3}} + a b c \left (\frac {a}{b}\right )^{\frac {1}{3}}\right )} \arctan \left (\frac {\sqrt {3} {\left (2 \, x - \left (\frac {a}{b}\right )^{\frac {1}{3}}\right )}}{3 \, \left (\frac {a}{b}\right )^{\frac {1}{3}}}\right )}{3 \, a b^{2}} - \frac {{\left (2 \, a e \left (\frac {a}{b}\right )^{\frac {2}{3}} + a d \left (\frac {a}{b}\right )^{\frac {1}{3}} - a c\right )} \log \left (x^{2} - x \left (\frac {a}{b}\right )^{\frac {1}{3}} + \left (\frac {a}{b}\right )^{\frac {2}{3}}\right )}{6 \, b^{2} \left (\frac {a}{b}\right )^{\frac {2}{3}}} - \frac {{\left (a e \left (\frac {a}{b}\right )^{\frac {2}{3}} - a d \left (\frac {a}{b}\right )^{\frac {1}{3}} + a c\right )} \log \left (x + \left (\frac {a}{b}\right )^{\frac {1}{3}}\right )}{3 \, b^{2} \left (\frac {a}{b}\right )^{\frac {2}{3}}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3*(e*x^2+d*x+c)/(b*x^3+a),x, algorithm="maxima")

[Out]

1/6*(2*e*x^3 + 3*d*x^2 + 6*c*x)/b - 1/3*sqrt(3)*(a*b*d*(a/b)^(2/3) + a*b*c*(a/b)^(1/3))*arctan(1/3*sqrt(3)*(2*
x - (a/b)^(1/3))/(a/b)^(1/3))/(a*b^2) - 1/6*(2*a*e*(a/b)^(2/3) + a*d*(a/b)^(1/3) - a*c)*log(x^2 - x*(a/b)^(1/3
) + (a/b)^(2/3))/(b^2*(a/b)^(2/3)) - 1/3*(a*e*(a/b)^(2/3) - a*d*(a/b)^(1/3) + a*c)*log(x + (a/b)^(1/3))/(b^2*(
a/b)^(2/3))

________________________________________________________________________________________

mupad [B]  time = 5.07, size = 319, normalized size = 1.56 \[ \left (\sum _{k=1}^3\ln \left (\mathrm {root}\left (27\,b^6\,z^3+27\,a\,b^4\,e\,z^2+9\,a\,b^3\,c\,d\,z+9\,a^2\,b^2\,e^2\,z+3\,a^2\,b\,c\,d\,e+a\,b^2\,c^3+a^3\,e^3-a^2\,b\,d^3,z,k\right )\,\left (6\,a^2\,e+\mathrm {root}\left (27\,b^6\,z^3+27\,a\,b^4\,e\,z^2+9\,a\,b^3\,c\,d\,z+9\,a^2\,b^2\,e^2\,z+3\,a^2\,b\,c\,d\,e+a\,b^2\,c^3+a^3\,e^3-a^2\,b\,d^3,z,k\right )\,a\,b^2\,9-3\,a\,b\,c\,x\right )+\frac {a^3\,e^2+b\,c\,d\,a^2}{b^2}+\frac {x\,\left (a^2\,d^2-a^2\,c\,e\right )}{b}\right )\,\mathrm {root}\left (27\,b^6\,z^3+27\,a\,b^4\,e\,z^2+9\,a\,b^3\,c\,d\,z+9\,a^2\,b^2\,e^2\,z+3\,a^2\,b\,c\,d\,e+a\,b^2\,c^3+a^3\,e^3-a^2\,b\,d^3,z,k\right )\right )+\frac {d\,x^2}{2\,b}+\frac {e\,x^3}{3\,b}+\frac {c\,x}{b} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((x^3*(c + d*x + e*x^2))/(a + b*x^3),x)

[Out]

symsum(log(root(27*b^6*z^3 + 27*a*b^4*e*z^2 + 9*a*b^3*c*d*z + 9*a^2*b^2*e^2*z + 3*a^2*b*c*d*e + a*b^2*c^3 + a^
3*e^3 - a^2*b*d^3, z, k)*(6*a^2*e + 9*root(27*b^6*z^3 + 27*a*b^4*e*z^2 + 9*a*b^3*c*d*z + 9*a^2*b^2*e^2*z + 3*a
^2*b*c*d*e + a*b^2*c^3 + a^3*e^3 - a^2*b*d^3, z, k)*a*b^2 - 3*a*b*c*x) + (a^3*e^2 + a^2*b*c*d)/b^2 + (x*(a^2*d
^2 - a^2*c*e))/b)*root(27*b^6*z^3 + 27*a*b^4*e*z^2 + 9*a*b^3*c*d*z + 9*a^2*b^2*e^2*z + 3*a^2*b*c*d*e + a*b^2*c
^3 + a^3*e^3 - a^2*b*d^3, z, k), k, 1, 3) + (d*x^2)/(2*b) + (e*x^3)/(3*b) + (c*x)/b

________________________________________________________________________________________

sympy [A]  time = 1.64, size = 178, normalized size = 0.87 \[ \operatorname {RootSum} {\left (27 t^{3} b^{6} + 27 t^{2} a b^{4} e + t \left (9 a^{2} b^{2} e^{2} + 9 a b^{3} c d\right ) + a^{3} e^{3} + 3 a^{2} b c d e - a^{2} b d^{3} + a b^{2} c^{3}, \left (t \mapsto t \log {\left (x + \frac {9 t^{2} b^{4} d + 6 t a b^{2} d e - 3 t b^{3} c^{2} + a^{2} d e^{2} - a b c^{2} e + 2 a b c d^{2}}{a b d^{3} + b^{2} c^{3}} \right )} \right )\right )} + \frac {c x}{b} + \frac {d x^{2}}{2 b} + \frac {e x^{3}}{3 b} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x**3*(e*x**2+d*x+c)/(b*x**3+a),x)

[Out]

RootSum(27*_t**3*b**6 + 27*_t**2*a*b**4*e + _t*(9*a**2*b**2*e**2 + 9*a*b**3*c*d) + a**3*e**3 + 3*a**2*b*c*d*e
- a**2*b*d**3 + a*b**2*c**3, Lambda(_t, _t*log(x + (9*_t**2*b**4*d + 6*_t*a*b**2*d*e - 3*_t*b**3*c**2 + a**2*d
*e**2 - a*b*c**2*e + 2*a*b*c*d**2)/(a*b*d**3 + b**2*c**3)))) + c*x/b + d*x**2/(2*b) + e*x**3/(3*b)

________________________________________________________________________________________